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Abstract 

Physical properties of the extracellular matrix (ECM) are known to regulate cellular processes ranging from spreading to 
differentiation, with alterations in cell phenotype closely associated with changes in physical properties of cells themselves. 
When plated on substrates of varying stiffness, fibroblasts have been shown to exhibit stiffness matching property, wherein 
cell cortical stiffness increases in proportion to substrate stiffness up to 5 kPa, and subsequently saturates. Similar 
mechanoadaptation responses have also been observed in other cell types. Trypsin de-adhesion represents a simple 
experimental framework for probing the contractile mechanics of adherent cells, with de-adhesion timescales shown to 
scale inversely with cortical stiffness values. In this study, we combine experiments and computation in deciphering the 
influence of substrate properties in regulating de-adhesion dynamics of adherent cells. We first show that NIH 3T3 
fibroblasts cultured on collagen-coated polyacrylamide hydrogels de-adhere faster on stiffer substrates. Using a simple 
computational model, we qualitatively show how substrate stiffness and cell-substrate bond breakage rate collectively 
influence de-adhesion timescales, and also obtain analytical expressions of de-adhesion timescales in certain regimes of the 
parameter space. Finally, by comparing stiffness-dependent experimental and computational de-adhesion responses, we 
show that faster de-adhesion on stiffer substrates arises due to force-dependent breakage of cell-matrix adhesions. In 
addition to illustrating the utility of employing trypsin de-adhesion as a biophysical tool for probing mechanoadaptation, 
our computational results highlight the collective interplay of substrate properties and bond breakage rate in setting de- 
adhesion timescales. 
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Introduction 

The extracellular matrix (ECM) which serves as a scaffold for 
maintaining the integrity of various tissues, is known to encode a 
diverse range of physical cues, including stiffness, topography, 
geometry, ligand spacing and dimensionality [1]. Of these, ECM 
stiffness has emerged as an important factor and has been shown 
to modulate a range of cellular processes including spreading [2] , 
motility [3], differentiation [4] and cancer invasion [5]. Such 
responses to ECM features require a close coupling between cell- 
matrix adhesions (also called focal adhesions) and the contractile 
acto-myosin cytoskeleton, leading to active reorganisation of the 
cytoskeleton and the adhesions [6,7]. Changes in cellular processes 
are closely tied to changes in physical properties of the cells, as 
evidenced by changes in cell cortical stiffness and traction forces 
exerted by a range of different cell types across substrates of 
varying stiffness [8]. 

Trypsin de-adhesion represents a simple method for probing the 
biophysical properties of adherent cells [9,10]. In this assay, upon 



incubation with warmed trypsin, cells round up driven by rapid 
severing of cell-matrix adhesions. The retraction process obeys 
sigmoidal kinetics with time constants that track cortical stiffness 
values. Further, de-adhesion time constants are sensitive to cellular 
contractility, with contractile activation leading to faster de- 
adhesion and contractile inhibition leading to delayed de- 
adhesion. Consequently, the de-adhesion assay has been used for 
studying modulation of cellular contractility by various features of 
the ECM. In breast cancer cells, increase in ECM density has been 
shown to increase protease-mediated ECM degradation in a 
contractility-dependent manner [11,12]. In addition to tracking 
cell rounding during de-adhesion, the pattern of cell movement 
during de-adhesion, i.e., translation and/or rotation, may differ 
from cell to cell, and was shown to depend on the spatial 
anisotropy in cell contractility, bond distribution and bond 
strength. Specifically, while asymmetry in bond strength and/or 
bond distribution was shown to cause cell translation, cell rotation 
required spatial asymmetry in both bond distribution and 
contractility [13]. These results may help us in understanding 
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the collective influence of contractility, bond distribution and bond 
strength in modulating random versus persistent cell motility. 
Taken together, these results illustrate the effectiveness of the 
trypsin de-adhesion assay for probing cell mechanics. 

Significant insight into the contributions of adhesion and 
contractility to various cellular processes has been obtained via a 
diverse range of theoretical and computational studies [14—16]. 
Several of these studies have tried to understand the role of 
substrate properties on cellular responses including spreading and 
motility [17-20]. Some of these studies include traction force 
localisation in cells [21], the role of substrate stiffness on stress fiber 
alignment [22], and the role of substrate thickness, stiffness and 
geometry of adhesion patches on traction forces generated by 
adherent cells [23,24]. Though a cell is a highly heterogeneous 
and dynamic entity, it is not uncommon in the literature to have a 
simple, linear, isotropic, elastic/visco-elastic description of the cell 
[25—28]. Using a similar description, Mofrad and co-workers 
demonstrated the viscoelastic behavior of NIH 3T3 cells during 
cell detachment [29,30]. To characterise focal adhesion dynamics, 
Schwarz and Erdmann had proposed a formulation where bond 
lifetime was dependent on cluster size, rebinding rate, and the 
force exerted at the adhesion [31]. In another study, Schwarz and 
co-workers used a two spring model to describe the contributions 
of cell-matrix adhesions to cellular activities, showing the influence 
of substrate stiffness on cell force built up [32]. Collectively, these 
studies indicate that in modelling cell-substrate interactions, cells 
and their underlying substrates are often modelled as elastic or 
viscoelastic solids, with simple description for capturing adhesion 
dynamics. 

While de-adhesion experiments have been successful in 
demonstrating contractile modulation by ECM density, it remains 
unclear if the same assay can also be used for probing the cellular 
mechanoadaptation response on surfaces of varying ECM stiffness. 
In this paper, we present a computational framework for probing 
the collective influence of cell and substrate properties, distribution 
of adhesions and cell contractility in regulating cell-matrix de- 
adhesion dynamics. Here, the cell is modeled as a pre-stressed 
viscoelastic solid adhering to an elastic substrate via cell-matrix 
adhesions at a prescribed density. De-adhesion is simulated by 
rate-dependent chopping of bonds. By studying the dynamics of 
retraction, we demonstrate how the time scales are regulated by 
substrate properties and bond breakage rate, and determine 
analytical expressions of the time constants for certain regimes of 
the parameters. While comparing our model predictions with 
experimental de-adhesion curves of NIH 3T3 fibroblasts on 
polyacrylamide substrates of varying stiffness, we observe that 
contrary to delayed de-adhesion on stiffer substrates predicted by 
our model, fibroblasts exhibit faster de-adhesion on stiffer 
substrates. Finally, we show that this anomaly can be resolved 
by introducing a force-dependent bond breakage rate in our 
model. Taken together, our results demonstrate the utility of 
employing trypsin de-adhesion as a biophysical tool for probing 
mechanoadaptation and illustrate the complex dependence of de- 
adhesion time constants on substrate properties. 

Materials and Methods 

Fabrication of polyacrylamide hydrogels 

Polyacrylamide hydrogels were fabricated as described else- 
where [33]. Briefly, glass coverslips were cleaned with 0.1 N 
NaOH, silanized with APTES for 5 mins (Sigma), rinsed with 
distilled water, and treated with 0.5% glutaraldehyde for 30 mins. 
Stocks of 40% acrylamide solutions and 2% bisacrylamide 
solutions (BioRad) were combined to obtain gels of 1.5 kPa, 



2.5 kPa and 40 kPa stiffness. Finally, gels were covalently cross- 
linked with 10 jig/ 'ml rat tail collagen I (Sigma) using the 
photoactivable bi-functional crosslinker sulfo-SANPAH (Pierce). 

De-adhesion experiments 

NIH 3T3 fibroblasts were maintained in T75 flasks in DMEM 
(Himedia) supplemented with 10% FBS (Invitrogen) at 37°C with 
5% C02. For experiments, fibroblasts were cultured on 
polyacrylamide hydrogels for a period of 24 hours. Quantification 
of cell spreading and cell shape was performed for at least 100 cells 
per condition across two independent experiments. For perform- 
ing de-adhesion experiments, cells were washed with phosphate 
buffered saline (PBS, Himedia), incubated with 0.25% warm 
trypsin (Himedia) and immediately imaged in phase contrast using 
an inverted microscope (Olympus 1X71) at 10 x magnification. 
De-adhesion experiments were performed on a temperature 
controlled stage with temperature set at 37°C. Cells were imaged 
till they became rounded but remained attached to the substrate. 
Movies were processed manually in Image J (NIH) to obtain 
experimental de-adhesion curves (i.e., plots of normalized change 
in area versus time), which were fitted with Boltzmann curve to 
determine the de-adhesion time constants T\/2 and 7*2 . De- 
adhesion time constants were obtained for 15 cells per condition 
across two independent experiments. 

Computational Model 

To understand the role of substrate stiffness on experimentally 
observed de-adhesion pattern, we developed a one-dimensional 
(1-D) model, with the variables being stiffness of cell and substrate, 
cell contractility, number and distribution of cell-substrate bonds, 
and rate of bond breakage. Unlike complex cellular processes such 
as cell spreading or cell migration which occur in the timescale of 
tens of minutes to hours [1,34,35], and involve complex 
cytoskeletal reorganisation, de-adhesion involves rounding up of 
a 50-100 jim. size cell in a matter of seconds, where perhaps, 
cytoskeletal dynamics are less important. This has motivated us to 
have a simple viscoelastic description for the cell, as done in many 
other studies [26,29,30,36-39]. Specifically, the cell was modeled 
as a 1-D viscoelastic object connected to the elastic substrate by a 
continuous distribution of bonds [32]. Various previous studies 
suggest that effects of cell contractility can be reasonably 
represented using pre-stress/pre-strain-based approaches [21,40— 
42]. For the present study, we have adopted a similar method to 
incorporate cell contractility. In order to mimic contractility in our 
model, the cell, initially in a relaxed state, was brought to a 
stretched state by the application of a specified displacement field 
and attached to the substrate below with the help of cell-substrate 
bonds (see Fig. 1). Once the cell is connected to the substrate, it 
starts applying tractions on the substrate, and ultimately reaches 
mechanical equilibrium (tensional homeostasis) with the substrate 
resisting the internal pre-stress/ contractile forces. Note that, in the 
absence of cell-substrate connections, no pre-stress buildup is 
possible. 

For simulating adhesion and subsequent de-adhesion, the 
stretched configuration of the cell (at time t = 0) was generated 
by imposing an initial displacement field (woO 7 )) to a cell of rest 
length Z/ re st [21], and connecting it to the substrate through cell- 
substrate bonds (Fig. 1(a)). The equilibrated length of the cell is Lq. 
Thus, at time t = 0, the cell is in a completely spread out state and 
in mechanical equilibrium with the substrate, the length co- 
ordinate being y. For this configuration, the stress along the cell 

can be written as o(y)=E c — where u c (y) is the cell 

displacement at any point along the length of the cell. Due to 
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Figure 1. Schematic representation of the model. Cell is modeled as a viscoelastic solid (represented by springs and dashpots), and is physically 
coupled to the underlying elastic substrate through focal adhesions drawn as springs. The number of springs depend on the density of adhesions, (a) 
Contractility is imposed by stretching the cell and then coupling it to the underlying substrate through cell substrate bonds, till mechanical 
equilibrium is attained. Inward directed arrows inside the cell indicate the cell pre-stress. (b) De-adhesion was simulated by breaking bonds in a time- 
dependent manner till the cell reached the final rounded configuration (time -*oo) where all the adhesions were broken. The spatial positions of the 
cell in the fully stretched (a) and completely relaxed (b) states are shown, (c) Relationship between different displacement components of the system 
at equilibrium is shown using a representative element P. 
doi:1 0.1 371 /journal. pone.01 0691 5.g001 



the interaction of the cell with the substrate, the stress is also 
related to the substrate displacement, u s (y). Thus, the force 
equilibrium equation along y direction can be written as 



My) | F _ Q 

dy bh c 



(i) 



where F is the force transmitted (per unit length) by cell substrate 
bonds, b is the width and h c is the height of the cell. To get a 
relation between the substrate displacement u s (y) and traction, we 
assume that the force f, which is transmitted by a single cell- 
substrate bond, is distributed uniformly over a circle of size df, the 
diameter of the bond (^10nm) [43,44]. Further, assuming the 
substrate with Young's modulus E s to be of much greater thickness 
as compared to the lateral extent of the cell, we use standard 
expressions to obtain a relation of the form / « E s diU s , where u s is 
the displacement of the center of the circle [45] . If the stiffness 
contributed from a single bond is E s df, then the overall 
contribution to the substrate stiffness from rib bonds is simply 
ribE s di. If Pq(j) represents the density of adhesions per length 



along the length of the cell and E s refers to the scaled effective 
stiffness of the substrate [43,44] given by the expression 
E s ribdi 

E s = — — — , the force equilibrium in Eq. 1 can be written as 



bh c 



do(y) 
dy 



-E s u s (y)p 0 (y) = 0, 



(2) 



After writing the substrate displacement u s {y) in terms of cell 
displacement as u s (y) = u^iy) — u c (y), and non-dimensionalizing 
all length dimensions with Lq, Eq. 2 becomes, 



d 2 u c (y) 
df 



+ n (u 0 (y) - u c (y))p(y) = 0, 



(3) 



where y = y/Lo, Po(y) — po(y)Lo and the dimensionless parameter 
E s Lq 

r\ = — - — . (In this paper, the convention followed for representing 
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symbols is that a quantity with dimension of length or time is 
denoted with corresponding symbol and any parameter with a () 
on top denote the corresponding dimensionless quantities of 
original equation). The displacement profile for the cell at / = 0 
can be obtained by solving Eq. 3, with stress free boundary 
conditions, (du c /dy)y = 0 ^=0. Previously, the bond density has 
been estimated to be of the order of 10 3 per square microns [46]. 
Since larger focal adhesions tend to be located at the cell periphery 
in many different cell types [47], significant rounding of the cell 
will occur only when these bonds are broken. Assuming a cell of 
10 jum radius and 1% adhesion belt at the periphery, one can 
estimate rib ~10 3 . For cells cultured on ligand-coated glass 
coverslips, taking cell stiffness E c = (l — 10) kPa and effective 
elasticity of ligand/glass complex to be 0.1 MPa [2], one can 
obtain estimates of r\ in the range of 10—1000. 

De-adhesion was initiated by time-dependent irreversible 
breakage of bonds at a rate r per unit time. Thus, upon chemical 
detachment of the bonds, the cell rounds up driven by the 
mechanical imbalance caused by the pre-stress in the cell. Finally, 
at time t=co, with all cell-substrate bonds broken, the length of 
the cell would be L CX) — L rest (fig. 1(b)). In our model, the stress 
along the cell at any time t is written as: 



<r(y,t) = E c 



du(y,t) , _ d 2 u(y,t) 



dy 



■ + G t 



dydt 



(4) 



where E c and G c are the elastic and viscous modulus of the cell, 
respectively. As in the earlier case, the stress in cell is also related to 
the substrate displacement. The equilibrium equation can be 
written as 



da(y,t) 

dy 



+ E s (u 0 (y) - u(y,t))p(y,t) = 0, 



(5) 



Taking the values of stress from Eq. 4, and re-writing substrate 
displacements in terms of cell displacements, Eq. 5 becomes as 



E c *^$^ + G c d ^ (6) 



dy 2 



d 2 ydt 



Non-dimensionalizing Eq. 6 using To(= G c /E c ) and Lq as 
metrics of time and length, respectively, we obtain the following 
non-dimensionalized version of Eq. 6: 



The time-dependent cell displacement u can be obtained by 
solving the above equation Eq. 7 using Eq. 8, subject to the 
following boundary conditions and initial condition: 



du(y,t) 



= 0, 



d y /j>=o,i 

u(y,Q) = u c (y). 



(9) 



In all our simulations, instead of separately providing values of 
E s , E c , and other variables, we vary one of the three parameters 
mentioned in table 1. 

Results 

De-adhesion of fibroblasts on hydrogels of varying 
stiffness 

Seminal work done by Janmey and co-workers demonstrated 
the mechanoadaptation response exhibited by fibroblasts cultured 
on substrates of varying stiffness [48] . Specifically, it was observed 
that fibroblasts match their cortical stiffness (E ce u) to that of 
substrate stiffness (E ge i) on substrates softer than 5 kPa, with 
cortical stiffness reaching saturation levels on substrates stiffer than 
5 kPa. To probe how substrate stiffness influences de-adhesion 
timescales, NIH 3T3 fibroblasts were cultured on 3 different 
substrates (1.5 kPa, 2.5 kPa and 40 kPa) [33]. Two of these 
substrates were chosen softer than the 5 kPa threshold below 
which E ce u = E ge i was observed by Solon et al. [48] ; the third 
substrate was chosen significantly stiffer than 5 kPa. Consistent 
with previous cell-on-gel studies [8,48], fibroblasts exhibited a 
stiffness-dependent spreading response with increased spreading 
on stiffer substrates (Fig. 2(b)). Cell shape was also found to change 
with cell stiffness, with cells exhibiting a more elongated 
morphology (i.e., low circularity) on stiffer substrates (Fig. 2(c)). 
While cells de-adhered slowest on the softest 1.5 kPa substrates, 
fastest de-adhesion was observed on the stiffest 40 kPa substrates 
(Fig. 3(a)). Fitting of experimental de-adhesion curves with Boltz- 
man equation allowed us to quantitatively assess the influence of 
substrate stiffness on de-adhesion time constants (Fig. 3(b)). While 
T\j2 dropped from 80 sec observed on 1.5 kPa gels to 30 sec on 
2.5 and 40 kPa gels, T2 dropped from 30 sec observed on 1.5 kPa 
gels to 15 sec on 2.5 kPa gels and 8 sec on 40 kPa gels (Fig. 3(c, 
d)). Taken together, these results indicate that similar to cell 
stiffness, de-adhesion timescales are sensitive to substrate stiffness, 
and may be used as metrics for probing mechanoadaptation. 



*™> + + n(m)-m,W P (y,D=o, (7) 

oy z d z ydt 



where t=t/xo, y=y/L 0 , p 0 (y) = p 0 (y)L 0 and r x = 



E s Lq 



The 



displacements are naturally non-dimensionalized with Lq. We 
have also assumed a simple rate dependence for bond breakage 
given by 



p(y,7) = p 0 (y)exp(-r 2 7). 



(8) 



This is a standard procedure to describe first order rate kinetics 
of bond breaking or any similar kinetic events. It simply means 

that, we have used an equation of the following form = — rp. variables (f2 = 0.5, t/o(p) = 0.2j> and Po(y)= I2(y — 1/2) 2 ) were 



Model predictions: Effect of substrate stiffness on 
de-adhesion timescales 

To quantitatively assess how different physical quantities 
influence de-adhesion dynamics, we tracked the length of the cell 
L(t) = L !X) +u(Lo,t) — u(0,t) as a function of time. As is done in 
experiments [9] , throughout the paper we present this length in a 

L Q -L(t) 

normalized fashion defined as C(t) = — — . To determine the 

M) _ L y 

influence of substrate stiffness on de-adhesion timescales, we varied 
the dimensionless parameter r\, which represents the ratio of 
substrate stiffness to cell stiffness, and is referred to as stiffness 
ratio. To assess the sensitivity of stiffness ratio in regulating de- 
adhesion dynamics, we systematically varied r\ over three orders 
of magnitude from 1 to 1000, and tracked the normalized de- 
adhesion profiles (Fig. 4(a)). In doing these calculations, all other 
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Table 1. Description of parameters used in the study. 





Variable 


Description 


TO 


G c 
E c 


r\ 


E S U 
E c 


ri 


n 0 



doi:1 0.1 371/journal.pone.01 0691 5.t001 



kept unchanged. Under these conditions, fastest de-adhesion was 
observed for r\ = 1 with a half saturation time constant %\/2 ~2to 
and slowest for r\ = 1000 with a half saturation time constant 
Ti/2~12to where %\j2 is the time when C{x\/2) =1/2. Though 
these results illustrate the influence of stiffness ratio on de-adhesion 
timescales, these predictions contradict our experimental findings, 
wherein faster de-adhesion was observed on stiffer substrates. 

Model predictions: Influence of adhesion strength and 
adhesion distribution on de-adhesion timescales 

De-adhesion is expected to be influenced by the strength of 
bonds formed between cells and their substrate. In addition, 
several experimental studies have highlighted the effect of 
substrate stiffness in influencing the formation and dynamics of 
focal adhesions [49,50]. To directly determine the influence of 
adhesion strength on de-adhesion profiles, we computed the de- 
adhesion profiles for different values of the non-dimensional bond 



(a) 



1.5 kPa 



2.5 kPa 



40 kPa 




1.5 2.5 40 
Stiffness (kPa) 



1.5 2.5 40 
Stiffness (kPa) 



Figure 2. Quantification of cell spreading area and cell shape 
on substrates of varying stiffness, (a) Representative phase 
contrast images of NIH 3T3 fibroblasts cultured on 1.5, 2.5 and 40 kPa 
polyacrylamide gels for 24 hours, (b) Quantification of cell spreading 
area. Cell spreading increased with increasing substrate stiffness. Red 
bars denote statistical significance (*p<0.05). (c) Quantification of cell 
shape (circularity) on substrates of varying stiffness. Cell elongation 
increased with increasing substrate stiffness. Red bars denote statistical 
significance (*p<0.05). 
doi:1 0.1 371 /journal. pone.01 0691 5.g002 



cutting rate ^£{0.5,1,2,5} while keeping all the other variables 
fixed. This rate is representative of the characteristic time required 
for breakage of adhesions, and, is therefore, indicative of adhesion 
strength of the bonds. As seen in Fig. 4(b), f2=0.5 led to a de- 
adhesion profile with a half saturation Ti/2~4to. A ten-fold 
increase in led to a significantly faster de-adhesion response with 
^1/2 ~ T 0- Subsequent increase in led to faster de-adhesion with 
saturation for r2>10. Moreover, the shape of the de-adhesion 
profile changed from a sigmoidal curve to that of a single 
exponential 1— exp( — 7) for high values of reflecting the 
internal relaxation time To = 1 of the cell. 

Though these calculations demonstrate the influence of bond 
breakage rate in setting de-adhesion time scales, in all the above 
calculations, the distribution of adhesions p^iy) remained the 
same. However experiments with a range of different cell types 
have demonstrated that cells form bigger adhesions in peripheral 
regions, and smaller sized adhesions towards the centre. To study 
the influence of distribution of adhesions, we compared the de- 
adhesion profiles arising from three different patterns of adhesion 
p 0 (y)e{\, 12(1/2- y) 2 , 80(1/2- yf} while keeping all the other 
variables constant (r\ = 10, =0.5, iio(y) = 0.2y). Further, the 
distributions of adhesions were chosen in a manner so as to ensure 
that the total number of adhesions remained constant. The 
resulting de-adhesion plots are shown in Fig. 4(c). In comparison 
to the first case where adhesion was uniform across the entire 
length of the cell, the other two distributions were non-uniform but 
symmetric about the cell centre with strongest adhesions at the two 
ends and no adhesion in the centre. As seen in Fig. 4(c), these two 
distributions yielded nearly the same de-adhesion profile suggest- 
ing that de-adhesion profile is independent of the adhesion 
distribution. An extreme limit of bond-distribution is when 
Po(y)= l/2(d(y) + d(y — 1)), where 3 is the Dirac-delta function; 
this is when the adhesions are concentrated only on the two 
extreme edges. In such a limit, we need to track only the edges of 
the cell, and the governing equation of motion will be reduced to 
an ordinary differential equation: 



U(t) + (1 + riiV 0 e ^) W (7) = ri N 0 (u 0 - u c )(l - 
u(0) = 0. 



rit ), 



(10) 



Here, u(i) is the cell displacement and No is the total number of 
bonds at equilibrium, respectively, at the end of the cell. All other 
symbols have their usual meaning. It can be seen from the dotted 
line in Fig. 4(c) that this limit places a right bound on the de- 
adhesion profiles. 

Analytical estimates of ti/ 2 and t/ 

Thus far, our simulations have allowed us to asses the individual 
influence of stiffness ratio (r\) and bond-breaking rate fa) on the 
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Figure 3. Trypsin induced de-adhesion of NIH 3T3 fibroblasts. 

NIH 3T3 fibroblasts were cultured on polyacrylamide hydrogels of 
stiffness 1.5, 2.5 and 40 kPa for 24 hours and then de-adhesion 
experiments were performed, (a) Upon treatment of trypsin, cells 
started rounding up. Sequence of time-lapse snaps of de-adhering cells 
on substrates of different stiffnesses are shown, (b) Plot of normalized 
area (defined as change in area at any time divided by net change in 
area during the entire de-adhesion process) as a function of time for 
cells de-adhering on different substrates. Fitting of experimental de- 
adhesion curves with Boltzman curve allowed us to determine the two 
time constants of de-adhesion- T i /2 and T 2 , respectively, (c) and (d) 
show the plot of T 1/2 and T 2 respectively, for de-adhering cells on 
three different substrate stiffnesses. Red bars denote statistical 
significance (*p<0.05). 
doi:1 0.1 371 /journal. pone.01 0691 5.g003 

de-adhesion response. Now we study their collective influence and 
identify different regimes where each of these parameters play a 
dominant role. From the computed de-adhesion profiles, two time 
constants -T1/2 and T/-were extracted, so that we could 
quantitatively study their dependence on r\ and rj. While T\n is 
the time needed to reach £ = 0.5, Tf is the time required to reach 
£ = 0.95 from £ = 0.5. Even though we see in previous section 
that de-adhesion time scales are sensitive to both r\ and V2, the 
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Figure 4. Influence of substrate properties, bond-breakage 
rate, and adhesion distribution on de-adhesion dynamics, (a) 

De-adhesion curves computed by varying stiffness ratio (n) from 1 to 
1000 keeping r 2l uoiy) and p 0 (y) fixed (see text), (b) De-adhesion curves 
obtained for bond breakage rate (^£{0.5,1,2,5}) keeping r\, u^iy), and 
p 0 (y) fixed, (c) Influence of distribution of adhesions on cellular 
retraction was assessed by comparing three different adhesion 
distributions (plotted in inset) while keeping all other variables 
unchanged. The dashed line corresponds to the case with adhesions 
at extreme ends of the cell, and gives the rightmost bound of the de- 
adhesion curves. 

doi:1 0.1 371/journal.pone.01 0691 5.g004 



extend to which T1/2 and Tf depend on r\ and remains unclear. 
To get some insights, T1/2 and Tf were plotted for increasing values 
of f2 for r\ — 100 (red points) and 0.5 (blue points) (Figs. 5(a) and 
(b)). Irrespective of the choice of n, increase in f2 led to faster de- 
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adhesion with fastest de-adhesion time scales of Ti/2= In 2 and 
Tf = In 10, respectively. Under certain limiting conditions for the 
parameters r\ and f2, we were able to obtain approximate closed 
form expressions for T\/2 and Tf. They are plotted as lines in the 
respective plots to demonstrate their accuracy of fitting. In the 

regime f\ > 10 and r2<l, the expression (red line, Fig. 5(a)) is 

seen to be a very good estimate of T\j2- The corresponding 

expression for Tf is given by and works well for < 1 (red 

line, Fig. 5(b)). For very low values of n, T1/2 is given by solution of 

the equation \+e~ l - 



1 —e~ rit = \ (Fig. 5(a), blue line) 

1 — f2 1—7*2 2 



and follows the computed solution very closely. For large f2, Tf 
converges to In 10 irrespective of r\ (Fig. 5(b), blue line). 

Model predictions: Influence of force dependent bond 
breakage on de-adhesion timescales 

Till now, we have addressed the influence of substrate 
properties and cell-substrate adhesions in regulating de-adhesion 
dynamics and shown that stiffness ratio and bond-breakage rate 
play important roles in setting de-adhesion time scales. Further, 
the analytical insights obtained from the model provide a 
quantitative understanding of how various physical parameters 
govern the de-adhesion profile. Though insightful, this model, 
however, failed to reproduce the experimental observations, 
perhaps due to its insensitivity to the magnitude of contractility 
(t/o), i.e., for any non-zero magnitude of contractility this model 
produces identical de-adhesion curves. This insensitivity is due to 
the linearity of the model in u(y,t) because of which, the 
magnitude Uq of contractility gets cancelled out from the 
numerator and denominator in the expression for normalised 
length C(t). As explained below, we address this shortcoming in 
the current model by coupling contractility and bond-breakage. 

Several theoretical and experimental studies in the past have 
demonstrated that bond lifetimes are sensitive to forces [51-54], 
To capture the effect of force dependence, it was assumed that 
bond breakage happens in a force dependent manner [55-57] 
given by the expression 



dp 
dt 



: -rexp[fp]p 



(11) 



where p is the bond density as in previous sections, and / 
represents the force acting on the bonds. All other variables have 
the same meaning as earlier. In order to understand the de- 
adhesion dynamics of cell for this force-dependent bond breakage 
conditions, simulations were done as earlier and corresponding de- 
adhesion profiles obtained for varying values of r\ (0.1-100) 
(Fig. 6). Interestingly, is seen to play an important role in 
determining the de-adhesion profiles. As the value of constant 
increases, the nature of de-adhesion profile also reverses. Figure 6 
(a)— (d) show that for lower values of de-adhesion happens faster 
on softer substrates. In contrast, when is very high, the nature of 
de-adhesion profile changes and fastest de-adhesion occurs for 
highest stiffness, which matches with the experimental results 
(Fig. 3). 

Model predictions: Influence of contractility on 
de-adhesion timescales 

In our mathematical model, contractility is introduced by 
applying a pre-stretch to the initially relaxed cell and then attaching 
the cell to the substrate below by cell-substrate bonds. This way, the 



cell is brought to a pre-stressed state before the beginning of de- 
adhesion process itself, with the applied displacement field, Uo(J>) 
being a measure of cell contractility, i.e., higher the value of Uo(y), 
greater will be the cell pre-stress. Thus, the effect of substrate 
stiffness on cell contractility is implicit in our model. The exact 
nature of this dependence can be seen by plotting the contractile 
force at equilibrium for different values of r\ and varying Uo(y). As 
seen in Figure 7(a), it is observed that for a given initial displacement 
field (uo(y)), the pre-stress (at / = 0) initially increases with substrate 
stiffness and then saturates. Moreover, higher contractility on higher 
stiffness substrates, or, higher contractility on identical stiffness 
substrates (effected experimentally by contractility-activating drugs), 
would predict faster de-adhesion. This is clearly seen in de- 
adhesion profiles plotted for three different contractility values 
(uo(y)e{0.05y, O.lj, 0.2y}) (Fig. 7(b)), where higher contractility 
drives faster de-adhesion. 

Discussion 

In this paper, we have combined experimental and computa- 
tional approaches to understand how substrate stiffness influences 
trypsin-induced de-adhesion of fibroblasts. First, by doing de- 
adhesion experiments on substrates of varying stiffness, we have 
shown that cells on stiffer substrates de-adhere faster, driven by 
higher baseline contractility, with blebbistatin treatment nullifying 
the stiffness-dependent effects. To directly probe the contributions 
of various factors to de-adhesion, a computational model was 
developed incorporating the effects of stiffness ratio, strength and 
distribution of adhesions, and the contractility. Analytical 
estimates obtained for certain regimes of stiffness ratio and bond 
breakage rate demonstrate the coupled dependence of de-adhesion 
timescales on the two parameters. Comparison of the computed 
stiffness-dependent de-adhesion profiles with the experimental de- 
adhesion curves suggest that faster de-adhesion on stiffer substrates 
is predicted by a force-dependent bond breakage rate. Collectively, 
our results highlight the direct influence of substrate stiffness in 
regulating de-adhesion time scales. 

Cell spreading, proliferation, differentiation and invasion have 
all been shown to depend on substrate stiffness, and are intimately 
tied to the mechanical state of the cell. Specifically, substrate 
stiffness has been shown to modulate the cortical stiffness of NIH 
3T3 fibroblasts in a profound manner, with cortical stiffness closely 
matching the substrate stiffness upto 5 kPa, and reaching a 
saturation value on stiffer substrates [48] . Similar mechanoadapta- 
tion responses has also been documented in cancer cells, and 
require both an intact actin cytoskeleton and focal adhesion 
proteins [8,58]. Since higher myosin II activity increases traction 
forces [59,60] and cortical stiffness measurements [8,61], we 
reasoned that trypsin-induced de-adhesion which is driven by cell 
contractility [9,10], can be used as an assay for probing cell 
mechanoadaptation. Consistent with our hypothesis, we show here 
that cells de-adhere faster in comparison to cells on softer 
substrates. The ease of conducting de-adhesion experiments, the 
simplicity of analysis makes it ideally suited for probing contractile 
mechanics and cellular mechanoadaptation responses in adherent 
cells. However, it remains to be seen if similar substrate stiffness- 
dependent de-adhesion responses are exhibited by other cell types. 

De-adhesion-induced rounding requires simultaneous breakage 
of cell-matrix adhesions and cell retraction. Consequently, de- 
adhesion timescales depend both on the initial state of the 
cystokeleton and the rate of bond breakage. This has been 
demonstrated in previous experiments where contractile stimula- 
tion via drugs including LPA and nocodazole, or faster dis- 
assembly of adhesions effected by using a higher concentration of 
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Figure 5. Analytical estimates of de-adhesion time scales. %\/ 2 and t/ are plotted as functions of r 2 for two values of r\\ 100 (red points) and 
0.5 (blue points). The corresponding colour lines are fits to various analytical expressions described in the text. 
doi:1 0.1 371 /journal.pone.01 0691 5.g005 



trypsin, both led to faster de-adhesion [9,10]. Such a coupled 
dependency of T\j2 on the two factors can be appreciated from the 
lnv\ 

analytical expression %\n — . For a constant bond cutting rate 

?2 

(i.e., constant r^), this equation predicts delayed de-adhesion of 
cells on stiffer substrates (i.e., increase in the value of r\). Though 
this prediction directly contradicts our experimental observations 
where cells on stiffer substrates de-adhere faster compared to those 
on softer substrates, the utility of the analytical expressions of de- 
adhesion timescales obtained for force-independent bond breakage 
lie in demonstrating the dependence of de-adhesion timescales not 
on stiffness directly but on logarithm of stiffness. This is why 
increase of stiffness from 1.5 kPa to 40 kPa reduces %\j2 by 50% 



only. At high bond cutting rates, the de-adhesion timescales are 
relatively insensitive to substrate properties. Together, these results 
suggest that the combination of substrate properties and bond 
severing rates determine de-adhesion timescales. 

Several experimental studies have demonstrated that both cell 
spreading and traction forces increase with substrate stiffness. This 
is made possible by the formation of larger adhesions, which are 
also more stable and experience larger forces. Additionally, 
experimental and theoretical studies on dissociation of bonds 
under force have demonstrated that dissociation rate of a bond 
increases exponentially under the presence of a force. This would 
suggest that in de-adhesion experiments, disruption of a single 
focal adhesion leads to enhanced load on the other bonds, thereby 
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Figure 6. Influence of force-dependent bond breakage on de-adhesion dynamics. Force-dependent bond breakage was simulated by 
varying the parameter, /?. Figure (a) - (d) show the de-adhesion profile for different values of p for varying stiffness ratio (0.1-100). While delayed de- 
adhesion on stiffer substrates was observed for low values of ^ (/? = 0,100), faster de-adhesion was observed on stiffer substrates at higher values of ft 

(£ = 5000). 

doi:1 0.1 371 /journal.pone.01 0691 5.g006 



PLOS ONE | www.plosone.org 



8 



September 2014 | Volume 9 | Issue 9 | el 0691 5 



Substrate Stiffness-Dependent De-Adhesion Dynamics 




Figure 7. Influence of cell contractility on de-adhesion dynamics. The effect of cell contractility on de-adhesion dynamics was studied for 
three contractility values: wo(y)e{0.05j>,0.1J>,0.2j>} at the cell end, while keeping other parameters fixed, (a) Plot of contractile force at equilibrium for 
varying stiffness ratios for different contractility values, (b) Influence of cell contractility on de-adhesion dynamics, with other parameters kept 
constant. On a substrate of given stiffness, higher contractility led to faster de-adhesion. 
doi:1 0.1 371 /journal.pone.01 0691 5.g007 



reducing their bond lifetimes. Force-dependent bond breakage has 
been incorporated in our simulations through the parameter ft. 
Similar to the force-independent bond breakage results, delayed 
de-adhesion on stiffer substrates was observed for low values of 
( = 100 and less). This trend started to get reversed for values of 
in the range of 1000, and, for high values of ft ( = 5000), similar to 
experimental observations, faster de-adhesion was observed on 
stiffer substrates. The simulation results thus suggest that the 
nature of stiffness-dependent de-adhesion response depends on the 
magnitude of contractile forces exerted by cells. It is conceivable 
that a delayed stiffness-dependent de-adhesion response may be 
observed in cells which are less spread and possess low baseline 
contractility. 

As is widely known, intracellular contractility is one of the key 
parameters which influences process of adhesion, as well as, the 
process of de-adhesion as studied here. Motivated by similar 
studies in the past [21,40,41], we have adopted a pre-stress-based 
approach to model cell contractility. In conjunction with substrate 
stiffness, the initial displacement field Uo(y) sets the magnitude of 
cell contractility, thereby accounting for the effect of substrate 
stiffness on cell contractility (Fig. 7). Consistent with previous 
studies [62], our model suggests a direct dependence of cell pre- 
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